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Abstract. Blazars with redshifts z < 0.1 are likely can- 
didates for detection at energies in the range 300 GeV - 
50 TeV with Cerenkov telescopes and scintillator arrays. 
We present 7-ray flux predictions for a sample of 15 nearby 
flat-spectrum radio sources fitting the proton blazar model 
of Mannheim (1993a) to their observed broad-band spec- 
tral energy distributions. At high energies, we use fluxes 
or flux limits measured by ROSAT, CGRO and the Whip- 
ple Observatory to constrain their spectra. We take into 
account absorption of the 7-rays by pair production with 
low energy photons of the diffuse infrared-to-optical pho- 
ton background produced by galaxies (cosmic absorption) 
and with low energy synchrotron photons of the blazar 
radiation field (internal absorption). Typically, the theo- 
retical spectra decrease much faster above TeV (photon 
index s w 3) than between GeV and TeV (s « 2) owing 
to internal absorption. 

The predicted fluxes are confronted with flux limits in 
the 20-50 TeV energy range obtained by the High Energy 
Gamma Ray Astronomy (HEGRA) experiment. Without 
cosmic absorption, the fluxes are about equal to the cur- 
rent sensitivity of HEGRA. Improved 7/hadron separa- 
tion techniques could render a detection by HEGRA pos- 
sible, if cosmic absorption by the far-infrared background 
at wavelengths ^100 fim is not exceedingly strong. 

Key words: radiation mechanisms: non-thermal - instru- 
mentation: detectors - BL Lacertae objects: general - cos- 
mology: diffuse radiation - gamma rays: theory 



1. Introduction 

One of the most remarkable results obtained with the 
spark chamber experiment EGRET onboard the Compton 



Gamma Ray Observatory is the discovery of a large num- 
ber of 7- ray p oint sources at hi gh ga lactic latitudes (Fich- 
tel et al. |1994| , Montigny et al. |l995p . Most of the sources 
have been identified as blazars, which comprise an impor- 
tant subclass of active galactic nuclei (AGN). Blazars are 
characterized by highly variable, featureless and polarized 
continuum emission. In the radio-to-soft X-ray range, the 
emission is commonly attributed to synchrotron emission 
from a collimated plasma jet (Blandford & Konigl 1979) 



presumably emerging from the rotating magnetosphere of 
an accreting superm assive black hole of mass M ~ 10 s - 
10 10 M o (Camenzind |l990j) . Due to relativistic bulk mo- 
tion, the radiation pattern of the jet is sharply peaked 
about its axis. In blazars, the jet axis is aligned closely 
with the line-of-sight, leading to Doppler-boosted emis- 
sion and superluminal motion of radio knots. 
Recent multifrequency observations of blazars (Maraschi 
et al . 1994 , Falomo et al., 1995| , Valtaoja & Terasranta 



1995|) show that mm, optical, X-ray and 7-ray lightcurves 



all trace the same event during a 'flare'. This favors 
a synchrotron-self-Compton (SSC) or proton-initiated- 
cascade (PIC, Mannheim et al. 1991) origin of the 7-rays 
at evolving shock fronts. Both mechanisms predict corre- 
lated TeV emission (Zdziarski & Krolik 1994, Mannheim 



1996) 



At high photon energies, the 7-radiation from remote 
cosmic sources is severely attenuated by pair production 
via 7 + 7 —>■ e + + e~ on low energy background photons 
(Gould & Shreder [19661 ). Stecker et al. fll992j ) propose to 
infer the diffuse near-infrared background radiation, which 
is difficult to observe directly, by measuring the expected 
TeV cutoff in the spectra of blaza rs at z « 0.1 owing to 
cosmic absorption. Salamon et al. ( |l994|) argue that 7-ray 
measurements could also be used as a yardstick to de- 
termine the cosmic distance scale. MacMinn & Primack 
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( 1996 ) point out that the diffuse infrared-to- ultraviolet 
background is directly related to models of galaxy forma- 
tion. Measuring the diffuse background from the infrared 
to the ultraviolet bands, requires observations of b lazars 
at 7-ray energies from 10 GeV (Madau & Phinney 1996) 
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Fig. 1. The 7-ray horizons t 77 (_E,z) = 1 corresponding to 
the two different diffuse background models shown in Fig. 2. 
The horizons were calculated assuming Q = 1, q a = 0.5, 
Ho ~ 75 km s _1 Mpc -1 and a photon number density evolv- 
ing with redshift as n'de' = (1 + z) 3 nde (conserved number of 
photons) 



to 50 TeV. Probing the Universe with 7-rays will allow for 
independent tests of cosmological scenarios based on the 
well-understood physics of pair production. 
Figure 1 shows the 7-ray horizon marking the relation be- 
tween photon energy E and redshift z of a 7-ray source 
defined by t 11 {E,z) = 1. The solid curve is obtained by 
numerically integrating the optical depth r 77 as given in 
Dwek & Slavin ( 1994 ) using the present-day background 
spectrum computed by averaging various g alacti c evolu- 
tionary scenarios of MacMinn & Primack ( 1996 ) shown 
in Fig. 2. This spectrum remains well below current ex- 
perimental limits for the near-infrared background (Billcr 
et al. 



1995) 



Practically, knowledge of the turnover en- 
ergy E for a given redshift would yield the value of the 
diffuse background photon density above the resonant en- 
ergy e = 2m c c 2 /[E(l + z) 2 ]. Compared to our standard 
scenario for E(z) shown as the solid line in Fig. 1, the 
turnover energy increases if (i) the Hubble constant has 
a value greater than 75 km s _1 Mpc -1 , (ii) the diffuse 
background radiation has a density lower than the value 
obtained by averaging the MacMinn & Primack models 
or (iii) the evolution of the diffuse background density 
n'(e')de' oc (1 + z) k is more shallow than k = 3 already for 
moderate redshifts (keeping SI = 1 and A = fixed). A 
low background photon density would indicate late galaxy 
formation, and thus only a small amount of cold dark mat- 
ter in the models of MacMinn & Primack. Effects (i)-(iii) 
can reduce the optical depth r 77 (i?, z) by a maximum fac- 
tor of ~ 3 below the values adopted in the present work 
assuming the extremal values H = 100 km s _1 Mpc" 1 , a 
present-day infrared background equal to the integrated 



number counts of IRAS galaxies (de Zotti et al. 1995) and 
k = 0. 

By inspection of Fig. 1 it is obvious that sources at red- 
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Fig. 1. Solid line: the infrared-to-ultraviolet diffuse back- 
ground radiation field adopted in the present work. Dashed 
line: a diffuse background assuming that the 7-ray spectrum 
of Mrk421 cuts off at TeV due to cosmic a bsorpt ion. Trian- 
gles denote estimates by Madau & Phinney ( 199£| ) of the op- 
tial-to-ultraviolet diffuse background based on deep galaxy sur- 
veys 



shifts less than ~ 0.1 are required to pin down the 7-ray 
horizon in the TeV range. Several low-redshift EGRET 
blazars have been found, of which Mrk421 and Mrk501 
have also been d etecte d at TeV energies (Punch et al. 
1992 , Quinn et al. 1996 ), and many more potential sources 
with similar overall spectrum exist. To extract information 
on the intergalactic absorption from observed TeV spectra 
requires detailed knowledge of the spectra intrinsic to the 
sources. 

In this paper, we therefore investigate a sample of nearby 
flat-spectrum radio sources by fitting the proton blazar 
model of Mannheim (1993a) to their multifrequency spec- 
tra. In Sect. 2 we briefly outline the model and explain the 
fitting procedure. In Sect. 3, we introduce the sample and 
proceed in Sect. 4 comparing the predicted fluxes above 
TeV with the sensitivity of HEGRA. Finally, we summa- 
rize our results on the feasibility of mapping the 7-ray 
horizon in the TeV range. 

2. Outline of the proton-initiated cascade model 

Recent multifrequency observations of blazars seem to 
support an explanation of their 7-ray emission in terms 
of relativistic shocks propagating down an expanding jet 
(at least this is a powerful and predictive worki ng hyp oth- 
esis). Multiple shocks (Courant & Friedrichs 1948) can 
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form behind a magnetic nozzle (van Putten 1995). Parti- 
cle acceleration at relativistic and oblique shocks by dif- 
fusive and drift mechanisms has been studied by Kirk & 
Hea vens (p989| ), S chneid er & Kirk ( p89| ), K irk fc Schnei- 
der (p89| ), W ebb (p89| ), Begelman fc Kir k ( |i~99"o| ), Kriills 
1990| ), Kirk ( |199S ), Ellison et al. ( |1995[ ). In the inner- 



most regions of the jet, at the beginning of an outburst, 
the shocks seem to emit an optically thin mm-to-optical 
synchrotron spectrum (in some cases reaching up to soft 
X-rays) with an associated X-to-7-ray spectral component 
which can dominate the emitted power. Further down the 
jet, the shocks become transparent to radio synchrotron 
emission so that the initial outburst shows up delayed and 
broadened in radio lightcurves. Morphologically, the radio 
outbursts seem to be associated with the birth of new 
knots in VLBI maps. 

It is tantalizing to identify the high energy spectral compo- 
nent with the inevitable synchrotron-self-Compton emis- 
sion from the traveling shocks. However, this would require 
that the photon energy density exceeds the magnetic en- 
ergy density in sources where the 7-ray luminosity exceeds 
the optical luminosity, i.e. it ra d 3> ub, since 



"rad 

u B 



(1) 



The same conclusion u rac j 3> «b follows from assuming 
that TeV photons are singly scattered optical synchrotron 



hv/(m e c 2 ) 



3.41CT 14 5 ± 7 2 and 



photons. Since x opt 

= f£opt7 2 : we obtain B ± = 3 10 13 xl pt /x tcv G 
310" 3 



' . = f£opt7~, we obtain ts± = 6 W^"x^ pt /x te 

G corresponding to the comoving-frame magnetic 
energy density of ~ 4 10~ 7 7j -2 erg cm -3 . This value is neg- 
ligible compared to the comoving-frame radiation energy 
density ~ 1 erg s" 1 (Sect. 3. 2). However, the energy den- 
sity of the relativistic electrons u IC \ should certainly exceed 
the energy density u r ad of the photons they produce im- 
plying u 10 i 3> kb . This raises the problem of how it re i can 
be maintained at such a high level in a quasi-stationary 
situation. The commonly invoked mechanism of diffusive 
acceleration does not operate at M ro i 3> "b , since particles 
diffuse by pitch-angle scattering off magnetic field fluctu- 
ations. The pressure of these fluctuations is expected to 
be near the saturation value usb ~ ub, thus imposing the 
condition u ro i < ub- 

This conceptual problem with the SSC mechanism is re- 
moved if it is assumed that protons instead of electrons 
scatter off the synchrotron photons. In interstellar space, 
the energy in relativistic protons largely exceeds the en- 
ergy in electrons which could also be true in radio jets. The 
reasons for this are unclear. However, if the proton accel- 
eration mechanism can tap the Max well tail of a thermal 
distribution (Malkov & Volk 1995 ), it is plausible that 
the proton energy density will generally be larger than 
the electron energy density. Protons in a high radiation 
density environment, such as the shocks presumably trav- 
eling down radio jets, produce ample pions and pairs by 
inelastic scattering off the electronic synchrotron photons 



if they reach high energies ~ 10 s 11 GeV (Biermann & 
Strittmatter 1987). 

The pairs and pions initiate electromagnetic cascades 
(Mannheim et al. 1991) at the shocks. For a simple Bland- 
ford & Konigl (1979) type jet, the predicted spectra match 
the observed spectra from radio frequencies to 7-rays very 
well (Mannheim 1993a). Contrasting any other model pre- 
dicting 7-rays from blazars, the proton blazar model im- 
plies source spectra reaching > TeV photon energies gen- 
erally. Internal absorption of the 7-rays above TeV by the 
low energy electronic synchrotron photons leads to a steep- 
ening of the spectra by one power above TeV (s « 3). 
The PIC flux varies simultaneously with the flux of the 
infrared-to-optical target photons and with the proton 
maximum energy. 

Denoting the emission volume as V (assumed to be equal 
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Fig. 3. Sketch of jet geometry and typical synchrotron + 
PIC spectrum. Solid lines show the emission from various r 
in the range Vb < r < r max . Local electronic synchrotron 
spectra are characterized with increasing frequency by a syn- 
chrotron-self-absorption turnover, an optical thin S u oc v~ Y ^ 
radio/mm spectrum, a S v oc v~ x infrared spectrum and a 
steep near-infrared/optical/UV spectrum. The superposition 
of many such spectra along the expanding jet yields a flat ra- 
dio spectrum with S v oc const, and a S y oc v^ 1 spectrum 
above the break frequency v^, turning over steeply at v c . The 
corresponding PIC spectrum rises roughly as S u oc ^ _ ( - 5-0 ' 9 ) 
(depending on the pair creation optical depth r 77 ) and steep- 
ens in the MeV to GeV range where S v oc v" 1 is typical. Above 
~ TeV, the PIC spectra steepen further by one power and then 
cut-off at ultra high energies. Dashed lines indicate the spec- 
trum emitted when a newly formed shock (£ ~ 1) enters the 
region at r = Vh where the target density for proton cooling 
rises sharply (the synchrotron radiation at target frequencies 
becomes optically thin at this point). Further down the jet, the 
7-ray luminosity decreases, since the proton maximum energy 
required for £ ~ 1 cannot be maintained by the nonlinearly 
evolving shocks 
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for protons and electrons), the energy density ratio of pro- 
tons and electrons as i] — u p /u c (77 w 100 in the Milky 
Way at GeV energies) and the ratio of electron and pro- 
ton cooling times as £ = t c /t p , the PIC luminosity can be 
expressed in terms of the SSC luminosity Eq.(l) by 

Lpic ^ u p Vtp 1 = rju e Vt~ 1 ~ TjLsyJctp 1 = rj£L syn (2) 

The cooling time ratio £ depends strongly on the max- 
imum particle energies, since t coo \ oc 1/7 where 7 de- 
notes the Lorentz factor of the particle. Balancing cooling 
and expansion time scales for protons and electrons (so 
that £ = 1) yields extremely large proton Lorentz factors 
7 P ~ 10 8 ~ n . In the case of £ = 1, the PIC luminosity can 
exceed the SSC luminosity by the, possibly large, proton- 
to-electron ratio 77. A crucial requirement is, of course, that 
acceleration is faster than expansion and that the Larmor 
radius of the gyrating particles does not exceed the radius 
of curvature of the shocks (Bell 1978). The latter crite- 
rion hurts protons more than electrons, since the proton 
Larmor radius is much larger than the electron Larmor 
radius. If the coherent shock structure is destroyed during 
the nonlinear evolution of the shock propagating down the 
jet, proton scattering at the maximum energy becomes 
inefficient (scattering length ~ Larmor radius in Bohm 
diffusion). We therefore expect that £(7;) oc 7 P , m ax/7c,max 
is decreasing with time. Multiple shocks generate some 
time-average (£) which will be relevant for our sample of 
randomly observed blazars. 

We compute the cascade spectra numerically as described 
in Mannheim et al. (1991) and Mannheim (1993a), but 
taking into account additionally Bcthe-Heitler pair pro- 
duction and proton synchrotron radiation assuming a 
photo-pair luminosity L e ± ~ 0.52^ (strictly valid only for 
straight a = 1 power law photon fields as a target; Sikora 
et al. 1987 ). The optical depth of the jet with respect to 
pair creation is computed as follows: Since the typical an- 
gle of emission in the comoving frame is 9' ~ tt/2, the 
intersection of the line of sight with the jet corresponds 
to its transverse radius r (denoted as r± in Mannheim 



1993a). The target radiation (>far-infrared) becomes op- 



tically thin at r = where we expect the largest 7-ray 
luminosity to be produced (see caption of Fig. 3). The 
comoving-frame optical depth is then given by 



rl s V n°'77 r b = au B e' 1 ln[i/ c /^ b ] V 77 fb 



where 



ub 



7j/3j0(l + 7/A c /A F 



(3) 



(4) 



(cf. Eq.(24) in Blandford & Konigl |1979| ). The pair- 
creation cross section reaches at the reso- 
nant target photon energy e' = 2(m e c 2 ) 2 / E' for a 7-ray 
with energy E'. Using the Lorentz invariance of the optical 
depth, i.e. t 77 (.E) = t' (E 1 ), we obtain the optical depth 
in the observer's frame. Synthetic cascade spectra were 



produced with the following input parameters: redshift z, 
jet Lorentz factor y 3 , angle to the line of sight 8, proton- 
to-electron cooling rate ratio £, proton-to-electron energy 
density ratio 77, jet luminosity L- Jet (magnetic field + rcla- 
tivistic particles), synchrotron cutoff frequency in the co- 
moving frame v' c . While £ affects only the cascade part 
of the spectrum, 77 also affects the electronic synchrotron 
spectrum for fixed Lj ct , since it imparts the given energy 
between electrons and protons. In general, the synchrotron 
and cascade spectrum are both rather robust. The most 
important effect is the dramatic variation of the observed 
flux for variations with the boosting angle 8. Variations 
of the cosmological parameters f2, A and H a from their 
adopted values 1, and 75 km s _1 Mpc do not signifi- 
cantly affect the model fits, they do affect the luminosities 
in Tab. 2 accordingly. 

3. Multifrequency modeling of nearby blazars 

3.1. The sample 

We investigated the spectra of blazars from the list of Fich- 
tel et al. (1994) which have a redshift z < 0.1 or unkown, 
and which are visible from the northern hemisphere. The 
sources typically have a flat or inverted radio spectrum 
> 1 Jy. We used published multifrequency data, ROSAT 
data from the all-sky survey as well as from public point- 
ings, preliminary flux limits from Whipple observations 
and HEGRA scintillator and AIROBICC flux limits. At 
the time of writing, 5 out of a total of 15 sources were 
detected by EGRET, 2 = Mrk 421 a nd Mrk501 w ere de- 
tected by W hipple (Punch et al. |l992| , Kerrick et al. |l995b , 
Uumn et al. |1996D , Recently, the detection of Mrk421 has 



been confirmed independently by the HEGRA Ccrcnkov 
telescopes (Petry et al. 1995). Extending the sample to 
z — 0.2 would roughly double the number of sources. 
To estimate the effect of intergalactic absorption, we as- 
sumed z — 0.3 for those sources, where no redshift was 
available. At z < 0.3, the corresponding host galaxies are 
generally expected to be seen. We used the mean diffuse 
background density shown in Fig. 2 for the computation of 
the absorption. 

The data sets are neither simultaneous nor complete, and 
the intrinsic source variability, which shows up as a large 
dispersion of the fluxes for a given photon energy (e.g. 
Te V outb urst of Mrk421 by a factor of — 10; Kerrick et 
al. 1995b| , Macomb et al. 1995), brackets the theoretical 
spectrum. Since dust-rich quasars such as 3C273 (Lichti et 
al. 1995, Mannheim 1993b) are not present in our sample, 
additional absorption of the 7-rays by infrared photons 
from the host galaxy can be neglected. Predictions of TeV 
fluxes are given within the above uncertainties which may 
well add up to a factor of ~a few. Table 1 summarizes the 
high energy measurements of the sources comprising our 
sample. 



K. Mannheim et al.: Beacons at the Gamma Ray Horizon 



5 



-13.0 : /• 




6.0 -4.0 -2.0 0.0 2.0 4.0 6.0 8.0 10.0 12.0 14.0 -4.0 -2.0 0.0 2.0 4.0 6.0 8.0 10.0 12.0 
logJE/eV] log, [E/eV] 



4.0 -2.0 0.0 2.0 4.0 6.0 
log 10 [E/eV] 



8.0 10.0 12.0 



Fig. 4. Solid lines: proton blazar model fits. Dotted lines: same model fits with external (cosmic) absorption taken into account 
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Table 1. Summary of high-energy data available for northern blazars with z < 0.1 or unkown redshifts. The numbers in the 
columns below the experiments give the observed differential energy flux or flux limit log 10 [E dN/dE (erg cm - s -1 )] at the 
threshold energy of the respective experiment 



Source 




Redshift 


ROSAT" 


CGRO 6 


Whipple c 


HEGRA d 


Tibet 6 


0116+319 


(MS 01166+31) 


0.059 


-12.3 


< -10.9 


< -10.9 


< 


-11.0 


— 


0430+052 


(3C120.0) 


0.033 


-10.2 


< -10.7 


— 


< 


-10.5 


— 


0716+714 




— 


-11.6 


-10.5 


< -10.8 




- 


— 


0802+243 


(3C192.0) 


0.0599 


— 


< -10.5 


— 


< 


-10.8 


— 


1101+384 


(Mrk421) 


0.0308 


-10.0 


-10.6 


-10.8 


< 


-10.5 


< -10.6 












AIROBICC: 


< 


-11.3 




1214+381 


(MS 12143+38) 


0.062 


-12.7 


-11.2 


— 


< 


-10.9 


— 


1219+285 


(ON 231) 


0.102 


-11.8 


-10.6 


< -10.9 


< 


-10.9 


— 


1404+286 


(OQ 208) 


0.0797 


-12.8 


< -10.7 


— 


< 


-10.8 


— 


1514+004 


(PKS 1514+00) 


0.053 


-11.9 


< -10.8 




< 


-10.2 




1652+398 


(Mrk501) 


0.033 


-10.3 


< -10.5 


-11.1 


< 


-10.9 














AIROBICC: 


< 


-10.8 




1727+502 


(II Zw 077) 


0.055 


-10.8 


< -10.7 


< -10.8 


< 


-10.5 




1807+698 


(3C371.0) 


0.0512 


-12.0 


< -10.5 










2200+420 


(BLLac) 


0.069 


-10.9 


-10.3 


< -10.9 


< 


-10.6 




2201+044 


(PKS2201+04) 


0.028 


-11.9 


< -11.1 




< 


-10.3 




2209+236 








-10.7 




< 


-10.7 





PSPC differential flux at 1 keV adopting a = 0.5 (S v tx v' 01 ) 

EGRET differential flux at 100 MeV adopting a = 1; Fichtel et al. (1994), Thompson et al. 1995, Montigny et al. (1995) 

For BL Lac, we give a preliminary EGRET flux (Hartman 1995) 

Cerenkov telescope data at a threshold energy of 0.35 TeV adopting a = 2; Kerrick et al. (1995a) 
Scintillator array data at the threshold 40 TeV adopting an exponentially decreasing spectrum 
For Mrk 421 and Mrk 501 we also give th e AIR OBICC data for E t h ^ 25 TeV 

Data are from Kiihn ( [l994| ), Karle et al. ( [l995| ) and Westerhoff et al. (1995b) on behalf of the HEGRA collaboration 
Differential flux from air shower data at 10 TeV adopting an exponential spectrum; Amenomori et al. (1994) 



3.2. Summary of results 

The data used to obtain the model fits are non- 
simultaneous and there are large gaps in frequency cov- 
erage. Therefore, the parameter choices, and hence the 
spectra, are not free of ambiguities. The accuracy of the 
model predictions cannot be much better than to a factor 
of ~a few. Nevertheless, some important results can be 
inferred from the model fits: 

Gamma-to- optical luminosity ratios: Among the randomly 
sampled nearby blazars, there is no source with a strong 7- 
flare. The 7-ray luminosities are L 7 « r)S,L opt ss L opt with 
a proton-to-electron ratio (77) = 80 + 30 and (log 10 [£]) = 
— 1.8 ± 0.5. The absence of flares could be due to either 
a low duty cycle for £ « 1 states or, alternatively, the 
acceleration process could generally be less efficient in 
BL Lacertae objects than in radio quasars (most nearby 
blazars in our sample are BL Lacertae objects). This could 
be related to the magnetic field structure which is pref- 
erentially perpendicular to the jet axis in low-luminosity 
BL Lacs and parallel in high-luminosity quasars (Bridle & 
Perley 1984). 

TeV j-rays: We note that promising candidates for detec- 
tion at TeV energies are 0802+243, 0430+052, 1514+004 
and 1219+285 (in addition to Mrk421 and MrkSOl). The 
TeV photon flux predicted for 0716+714 accounting for 



cosmic absorption at z = 0.3 is similar to that for 
1219+285. However, we consider this redshift as a lower 
limit implying that the cosmic absorption is likely to be 
much stronger than assumed. 

Spectral slope between TeV and 50 TeV: Above TeV, 
the spectra continue as power laws, but with a steeper 
slope. The intrinsic source spectrum is given by I s = 
J (l — exp[— r 77i mt])/T 77 ,int corresponding to absorption 
in a homogeneous source, for which the photospheric 
size is energy-dependent. Hence, for T 77i j nt 3> 1 and 
since r 77 ,i nt cx E, the asymptotic spectrum is given by 
I s oc I /E. This is in marked contrast to the external 
absorption, which is given by an exponential law, i.e. 



,exp[ 



' 77.cxtJ 



Similarity of comoving-frame spectra: Somewhat surpris- 
ing, we find that the comoving-frame break frequencies 
are roughly constant with v' h ps 10 11 Hz. The comoving- 
frame cutoff frequencies also agree with a constant value of 
v' c w 3 10 14 Hz, which is theoretically expected from non- 



relativistic theory (Biermann & Strittmatter 1987). This 
could be a hint that the physics in the comoving frame 
of the relativistic inner jet is basically the same as in the 
(non-relativistic) hot spots of the jet hundreds of kilopar- 
secs away from the origin of the jet (Harris et al. 1994). 
Proton maximum energy: The average proton maximum 
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Fig. 5. Integral fluxes and cumulative HEGRA flux limit of the 13 blazars visible at La Palma (excluding 0716+714 and 
1807+698) from 10 GeV to 100 TeV. Dashed bold line: sum of fluxes with internal absorption. Solid bold line: sum of fluxes with 
internal+external absorption 



energy in the observer's frame inferred from the fits is 
~ 10 10 GeV. Detection of the resulting neutrino back- 
ground from blazars seems possible (Mannheim p.995 ). 
Jet Lorentz factor: The average jet Lorentz factor 
(7j)obs = 10 ± 4 is consistent with the value 7j = 7 in- 
ferred for low-luminosity radio sources based on the uni- 
fication of BL Lac s with Fanaroff-Rilcy type I galaxies 
(Urry & Padovani 1995| ). A freely expanding jet should 
have $ ss Tj -1 ss 6 deg which is consistent with the value 
from the fits. 

Variability time scales: From the model fits and Eq.(4) we 
obtain the typical numbers (u syn /u B ) — 0.01, (B) = 40 G, 

10. Using 

(6)(E>) = 



(rb) = 10 15 cm and the Doppler factor (i5) : 
Eq.(3), we obtain the 7-ray break energy (E) 
1 TeV where intrinsic absorption sets in (r^, 7 (i?') = 1). 
These numbers imply a minimum variability time scale 
(At ohs ) = (r b )(l + z)/((S)c) w 310 3 s in the observer's 
frame. 



4. Current flux limits 

In the TeV energy region, the search for point sources 
is only possible using earth-bound experimental setups 
sensitive to the secondary particles induced by cosmic 
ray primaries in the Earth's atmosphere. The main task 
is to separate the air showers induced by primary 7's 
from the much more abundant hadronic showers. Apart 
from detecting the charged secondary particle component 
using a matrix of scintillator counters, the photons of 
the Cerenkov light cone produced by relativistic electrons 
also give valuable complementary information about the 
shower and thus the incoming primary particle. Whereas 
up to now, scintillator arrays have provided only upper 
limits for the flux from 7-ray point sources, the employ- 
ment of imaging Cerenkov telescopes has been the most 
successful ground based observation technique of the last 
years: observation of ultra high energy 7-radiation from 
the Crab nebula, Mrk421 and Mrk501 are reported. As an 
example for experimental setups searching for TeV 7-rays, 
we quote the results obtained by the HEGRA detector ar- 
ray located on the Canary Island La Palma (28.8 N, 17.7 
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Table 2. Physical parameters for the model fits shown in Fig. 4 and predicted high-energy 7-ray fluxes. The last four columns 
give the predicted integral 7-ray fluxes above 1 and 40 TeV in units of 10" 1 cm" s" 1 accounting for internal absorption (index 
i) and (ii) internal + external absorption (index a) adopting the same differential spectra as noted in Tab.l 



Source 


7j 


6 
deg 


4> 

deg 


iogio 


[ hei. ' 
crg/s 


V 


1o Sio^J 


logio 


v c 
Hi 


r\\j> lev ) 


i* a ^> lev ) 


-Ti^> 4U le v ) 


-P a (,> 4UleV ) 


0116+319 


6 


6 


7 


44.8 


110 


-2.1 


14.8 


50 


30 


0.5 


310 -2 


0430+052 


8 


8 


5 


45.0 


30 


-1.3 


13.8 


600 


400 


0.2 


510~ 2 


0716+714 


16 


2 


1 


44.5 


70 


-1.3 


14.4 


900 


80 


20 


6IO- 9 


0802+243 


4 


3 


2 


44.0 


80 


-1.5 


14.1 


500 


500 


4 


210" 1 


1101+384 


20 


2 


5.7 


45.2 


50 


-2.2 


15.8 


500 


400 


0.6 


no -1 


1214+381 


5 


15 


5 


4D. t 


100 


-2.2 


1/19 
14. Z 


50 


30 


0.6 


210~ 2 


1219+285 


10 


7 


4 


45. z 


100 


-1.8 


1 A O 

14.0 


200 


80 


0.2 


710~ 4 


1404+286 


10 


3 


4 


44.7 


100 


-2.5 


13.7 


40 


20 


0.4 


710 -3 


1514+004 


7 


5 


3.5 


44.5 


100 


-2.0 


14.2 


200 


200 


0.6 


410" 2 


1652+398 


10 


2 


3 


43.9 


100 


-2.5 


16.5 


100 


100 


1.8 


610" 2 


1727+502 


10 


5 


10 


44.0 


30 


-1.1 


13.8 


80 


60 


0.08 


610" 3 


1807+698 


10 


9 


5.7 


44.9 


30 


-1.8 


14.1 


4 


3 


0.02 


1 10 -3 


2200+420 


10 


7 


10 


45.6 


100 


-1.9 


13.8 


50 


30 


0.2 


610~ 3 


2201+044 


10 


5 


4 


44.0 


90 


-2.2 


14.3 


70 


60 


0.2 


410~ 2 


2209+236 


10 


7 


5 


45.7 


100 


-0.8 


13.6 


400 


40 


0.7 


no- 10 


mean 


10 


6 


5 


44.8 


80 


-1.8 


14.4 










variance 


4 


3 


3 


0.6 


30 


0.5 


0.7 











W, 2200m a.s.L). As it comprises a scintillator array of 
224 counters, 17 Geiger towers for 7/hadron separation, 
and a 7 x 7 matrix of open Cerenkov counters (AIRO- 
BICC), HEGRA is able to equally exploit the information 
provided by the Cerenkov photons and the charged par- 
ticle component of air showers. The observation of single 
objects is furthermore possible by the use of 3 Cerenkov 
telescopes which form the first part of a system of 5 tele- 
scopes. The operation of AIROBICC is restricted to clear, 
moonless nights, but the solid angle acceptance of about 
1 sr allows to observe a large number of sources simulta- 
neously. As AIROBICC and scintillator arrays are sensi- 
tive to different physical properties of air showers, they 
have different energy thresholds: the scintillator array de- 
tects showers with more than about 40 TeV energy of 
the primary particle, the threshold for AIROBICC is at 
about 20 TeV. As a consequence of the lower Cerenkov 
light yield and larger fluctuations in proton showers, the 
threshold energy for proton induced showers exceeds the 
7-shower thresholds by about 10 TeV. The upper flux lim- 
its for 7-point sources from the AIROBICC array are of 
the order 10 _13 cm- 2 s" 1 with 7-shower energy thresholds 
of about 25 TeV. For the extragalactic sources Mrk421 
and Mrk501 AIROBICC provides upper flux limits of 
1.2410" 13 cm- 2 s- 1 and 3.72 10~ 13 cm" 2 s" 1 at threshold 
energies of 25.3 and 24.0 TeV respectively (Karle et al. 
1995). For the scintillator array, Kiihn (1994) gives an 
upper flux limit of 4.3 10~ 13 cm" 2 s _1 above 40 TeV for 
Mrk421. The results quoted above were obtained without 
any attempt to reduce the large background of hadron in- 
duced showers, which are about 10 5 times more numerous 
than 7-induced showers. 

Current data analysis applying cuts on the basis of 



7/hadron- separation techniques (Westerhoff et al. 1995a) 
yields a cumulative flux limit of 710~ 14 cm" 2 s" 1 (90% 
CL) above 50 TeV for the sum of the sources in our sam- 
ple excluding 0716+714 an d 1807 +689 with poor visibility 
at La Palma (Westerhoff 1995b). This is approximately 
equal to the sum of the unabsorbed theoretical fluxes de- 
rived from the model fits (Fig. 5, Tab. 2). Accounting for 
absorption, the theoretical flux drops by roughly an order 
of magnitude lower. It is difficult to assess the systemat- 
ical errors in the determination of the experimental flux 
sensitivity, but they can be estimated to be quite large. 
Since the theoretical errors are also difficult to control bet- 
ter than within factors of ~a few by the very nature of the 
variable blazars, a detection of blazars at ~50 TeV does 
not seem beyond reach. On the other hand, v ery st rong 
absorption, such as pro posed by Stecker et al. ( 1992j ) and 
Dwek & Slavin ( 1994 ) indicated as the dashed lines in 
Figs. 1 and 2, would further reduce the predicted flux 
above 50 TeV by orders of magnitude. In this case, de- 
tection of blazars at HEGRA scintillator energies would 
be truly impossible. 

Future analysis will include methods to increase the 
signal-to-background-ratio by fully exploiting the infor- 
mation of all components of the experimental setup and 
by employing new and improved methods of data analysis. 



5. Conclusions 

The continuum spectra of 15 nearby blazars have been 
modeled satisfactorily over ~ 20 orders of magnitude 
in frequency combining proton-initiated unsaturated syn- 
chrotron cascades and electron synchrotron radiation 
emitted by relativistic jets with a mean Lorentz factor 
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(7j) = 10 ± 4. The 7-ray spectra are complex, but a mean 
power law photon index s w 2 represents the spectra in the 
EGRET range rather well. The inferred physical param- 
eters are in line with a proton-to-electron energy density 
ratio (77) = 80 ± 30. The proton cooling rate is less than 
the electron cooling rate by a factor ~ 0.016t^oii imply- 
ing about equal infrared-to-optical and 7-ray luminosities 
from the same emission volume. 

Internal absorption of the 7-rays by the low energy syn- 
chrotron photons leads to a steepening of the spectra at 
- TcV above which s « 3 up to - 100 TcV. The TcV 
spectra of Mrk421 and Mrk501 are not significantly af- 
fected by cosmic absorption and should have s ~ 3.2. 
Cosmic absorption must be present in 0716+714 for which 
the predicted TeV emission without absorption greatly ex- 
ceeds the Whipple limit. Cosmic absorption becomes suf- 
ficiently strong for 0716+714 if its redshift is greater than 
0.3. 

Without cosmic absorption, the sum of the blazar fluxes is 
about equal to the current sensitivity of the HEGRA ar- 
ray. Taking into account cosmic absorption based on mod- 
els of galaxy formation by MacMinn & Primack (1996), 
the predicted flux falls an order of magnitude below the 
current flux limit. The absorption at ~ 50 TeV is not very 
sensitive to cosmological parameters (e.g., fl or A). A Hub- 
ble constant of 100 km s _1 Mpc -1 and late-epoch galaxy 
formation producing an infrared background equal to that 
produced by IRAS galaxies would increase our theoretical 
7-ray fluxes by a factor of ~ 3. A mild enhancement of 
the theoretical flux is expected from forward cascading of 
the pairs produced in the intergalactic medium by inverse- 
Compton scattering off background photons (Prothcroe & 
Stanev 1993, Aharonian et al. 1994, Plaga 1995), if the in- 
tergalactic magnetic fields are weak. 

A critical appraisal of the systematic errors in the deter- 
mination of the exact HEGRA energy threshold and flux 
level as well as the variable nature of the blazars leads us 
to conclude that detection of the nearest blazars with the 
HEGRA array seems feasible. A positive detection would 
disprove any claims of an exceedingly strong cosmic in- 
frared background (with profound implications on the era 
of galaxy formation) and would provide strong support 
for a baryon-induced cascade origin of the 7-rays from 
blazars. 
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Note added in proof: After submission of the manuscript, we 
became aware of a preprint by Puget et al. (1996, A&A, in 
press), claiming the detection of a diffuse FIR background us- 
ing COBE data. Within the uncertainties, their measured flux 



is equal to the far-infrared flux of the mean MacMinn & Pri- 
mack background spectrum used in our paper. 
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